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Using a new form of density functional theory for the ab initio description of electronic systems 
in contact with a dielectric environment, we present the first detailed study of the impact of a 
solvent on the surface chemistry of Cr2 03, the passivating layer of stainless steel alloys. Compared 
to vacuum, we predict that the presence of water has little impact on the adsorption of chloride 
ions to the oxygen-terminated surface but a dramatic effect on the binding of hydrogen to that 
surface. These results indicate that the dielectric screening properties of water are important to the 
passivating effects of the oxygen-terminated surface. 



1. INTRODUCTION 

Ab initio calculations have shed light on many physico- 
chemical questions, including chemical reactions in so- 
lution and chemical reactions at surfaces 

USB ai- 

though computational chemistry is now able to provide 
not only qualitative but also quantitative insights into 
surface chemistry |2j, so far these studies have been lim- 
ited to reactions in a vacuum even though experimentally 
these reactions often occur in a solvent environment. 

Each year, corrosion costs the United States $276 bil- 
lion 0, approximately 3.1% of the gross domestic prod- 
uct, and many approaches have been used to under- 
stand and model this complex process [E IS 0- High 
performance stainless steel alloys contain chromium to 
result in the formation of a Cr203 passivating surface 
layer, which provides corrosion resistance. Even with 
this layer, such alloys are susceptible to breakdown in 
acidic, chlorine-containing aqueous environments[8l |9j, 
the study of which demands simultaneous treatment of 
surfaces, reactants and the aqueous dielectric environ- 
ment. Direct experiments are difficult, and although ab 
initio calculations have a history of answering many such 
questions, to date they have not been able to address 
the effect of the solvent on surface reactions. Alavi et 
aL^(|j f° r instance, have studied the adsorption of HC1 
on single-crystal alpha- AI2O3 (0001) surfaces and calcu- 
lated adsorption energies as a function of surface coverage 
within density-functional theory, but all of their results 
are obtained in a vacuum environment. 

Previous studies of Cr 2 03 have been limited to pure 
surfaces in vacuum |ll|. an unrealistic environment for 
the study of corrosion. These studies indicated that 
the highest occupied and lowest unoccupied molecular 
orbitals are localized on the chromium ions, suggesting 
that the oxygen-terminated surface could provide a sta- 
ble barrier against acidic chlorinated environments if the 
surface oxygen layer could be prevented from reacting 
with species in the solution. The work below puts forth 
evidence to support the novel hypothesis that the dielec- 
tric screening effects associated with an aqueous envi- 
ronment actually prevent the formation of bonds with 
aqueous species such as protons, thereby rendering the 



oxygen-terminated surface virtually non-reactive. 

Models for calculating solvation energies from contin- 
uum dielectric theorv|12j have been applied to single 
molecules or activated complexes but not to molecules 
adsorbed on surfaces, perhaps because such methods are 
generally applied to molecules. Here, we introduce the 
first approach to ab initio calculations in a dielectric en- 
vironment which sits on a firm theoretical foundation. 
Below, we show that this new approach, which in a sim- 
ple approximation is related to that recently introduced 
by Fattebert and Gygi^jjJ, gives results in good agree- 
ment with currently accepted quantum chemical meth- 
ods and is well-suited to surfaces. Finally, we apply the 
approach to carry out the first ab initio study of the reac- 
tivity of hydrogen and chlorine on an oxygen-terminated 
Cr2O3(0001) surface in contact with a solution. 

2. THEORETICAL APPROACH 

Because of the large cells required in this study (well 
over one hundred atoms), the only practicable ab initio 
approach suited to describe the electrons of the surface 
and the reactants is density-functional theory, 11]. The 
periodicity of the surface makes periodic boundary con- 
ditions and thus the plane- wave pseudopotential method 
|l4| the most natural choice. The size of the surface cell 
and time-scales needed for proper thermodynamic aver- 
aging make a direct molecular dynamics treatment of the 
aqueous environment infeasible, thus raising the question 
of how to treat the solvent. This work takes the novel ap- 
proach of exploiting the existence of so-called "classical" 
density- functional theories^^, which can treat water rig- 
orously in terms of a simple thermodynamically averaged 
molecular density and an approximate functional|l6j]. We 
then introduce here for the first time the concept of a 
joint density-functional theory (JDFT) between the elec- 
trons in the surface and the molecules comprising the 
solvent. 

In a companion to this paper [l7^. we prove that the 
total thermodynamic free energy A of an electronic sys- 
tem (solute) with nuclei of charges Zj at locations Ri in 
equilibrium with a closed-shell liquid molecular environ- 
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ment (solvent) may be determined from the variational 
principle 

A = mm (A KS [n(r),{Zj,Rj}] (1) 

n(r),N(r) 

+A lq [N(r)} + U[n(r), N(r), {Z^Ri}]) , 

where, at the minimum, n(r) is the thermo- and 
quantum- mechanically averaged density of the elec- 
trons of the solute and N(r) is the likewise averaged 
molecular density of the solvent. AKs[n{r), {Zi, Ri}], 
Ai q [N(r)] and U[n(r), N(r),{Z I: R r }] appearing above 
are, respectively, the standard Kohn-Sham electron- 
density functional of the solute when in isolation, the 
classical density-functional for the molecular solvent 
when in isolation, and a new functional describing the 
coupling between the systems. The new functional 
U[n(r), N(r), {Zj , Ri}) is universal in the sense that it 
depends only upon the nature of the solvent and is inde- 
pendent of the nature of the solute. For completeness, we 
note that each of the three aforementioned functionals is 
implicitly also a function of the temperature T. 

In the present work, we make the further simplifica- 
tion of performing the minimization over N(r) in ifl}. 
resulting in the variational principle 

A = mmiAKsH^^Zr^Ri}} (2) 

n(r) 

+W[n(r),{Z I ,R I }]), 

where 

W[n(r),{Zi,Ri}] = min(A k [N(r)} 

N(r) 

+U[n(r),N(r),{Z I ,R I }]) 

is a universal functional dependent solely upon the iden- 
tity of the environment (and, implicitly, the tempera- 
ture). Note that, in principle, the theory at this stage 
is exact. Below we outline the approximations which 
we introduce because the exact form of the functionals 
AKs[ n i r ), {Zi, Ri}] and W[n(r),{Zi, Ri}} are unknown. 

3. COMPUTATIONAL DETAILS 

For treatment of the electrons in the chromium-oxide 
surface through the functional AKs[n(r), {Zj, Ri}], we 
apply the standard local spin-density approximation 
(LSD A) 01- The calculations themselves employ the 
total-energy plane-wave density-functional pseudopoten- 
tial approach fi3 with potentials of the Kleinman- 
Bylander form [l9( with p and d non-local corrections at 
a cutoff of 40 hartrees. Supercells with periodic bound- 
ary conditions in all three dimensions represent the sur- 
faces of isolated oxygen-terminated (OOOl)-oriented slabs 
of Cr 2 03 of thickness 13 A separated by 7.8 A of either 
vacuum or solvent. The in-plane boundary conditions 
suffice to describe 2x2 reconstructions and consist of four 
times the unit from Jll.] so as to allow sufficient isolation 
of solvent volumes excluded by chlorine adsorbed on the 



surface. The supercell contains a total of forty chromium 
and seventy-two oxygen atoms with chlorine or hydrogen 
added to the two surfaces in inversion-symmetric pairs, 
one on each side of the slab. Finally, we use a single 
k-point to sample the Brillioun zone of the surface slab. 
Ref. ^3 establishes that this choice of functional, pseu- 
dopotential, plane-wave cutoff, sampling density in the 
Brillioun zone and supercell gives a good description of 
the bulk and surface of Cr203. As in the aforementioned 
work, we employ the analytically continued functional 
approach [l4l l20j to minimize the Kohn-Sham energy 
with respect to the electronic degrees of freedom. Below, 
we relax all ionic configurations until the total energy is 
within 0.027 eV of the minimum and the maximum force 
in any direction is less than 0.3 eV/A. 

For the environment functional W[n(r), {Zi , Ri}] in 
(|2J), we take the interaction of the electron and nuclear 
charges of the system under study with a dielectric envi- 
ronment in which the dielectric constant is local in space 
and has a value dependent only upon the electron density 
at each point, e(r) = e(n(r)). In practice, we compute 
the total free energy of the system in contact with the en- 
vironment by finding the stationary point with respect to 
both the electrons and the mean electrostatic field (f>(r) 
of the functional 

A = A T xc[n^r), ni (r)] + Ay ps [n t (r),nj(r)] (3) 
+ J d 3 r L(r) U tot (r)-^^ 3 )(r-i^)j 

-^W)| 2 } 

where n^(r), n±(r) and ntot{r), respectively, are the up-, 
down- and total electron densities, ATxc[ n i( r ), n l( r )] is 
the Kohn-Sham single-particle kinetic plus exchange cor- 
relation energy within the local spin-density approxima- 
tion, AVps is the difference in the total pseudopotential 
energy from that expected from pure Coulomb interac- 
tions with point ions of valence charges Zi at locations 
Ri, and S^(r) is the three-dimensional Dirac-<5 function. 
Note that, although we do work directly with the Kohn- 
Sham orbitals, for brevity we have written the above in 
terms of the electron densities. We also note that at this 
level of approximation our joint density functional the- 
ory takes the same form as the approach of which 
introduced the form as a computational device without 
formal justification or proof. 

For the local dielectric function e(n(r)), we choose a 
specific form which varies smoothly from the dielectric 
constant of the bulk solvent e?, when the electron density 
n(r) is less than a critical value n c indicative of the in- 
terior of the solvent to the dielectric constant of vacuum 
e = 1 when n(r) > n c . Specifically, we take 

£( „ H1 + ^ ertc (M^>), (4) 

where the parameter a, to which the results are not very 
sensitive, controls the width of the transition from bulk 
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Figure 1. Comparison of predictions of a simple joint 
density-functional theory (vertical axis) with 
established quantum chemical values (horizontal axis) 
for acetamide, acetic acid, methanol, ethanol and 
methane (from left to right) from [2lj . 



to vacuum behavior. To determine the parameters a 
and n c for this simple model, we fit computed electro- 
static solvation energies in aqueous solution (e& = 80) 
to the accepted values from the quantum-chemical lit- 
erature for methane, ethanol, methanol, acetic acid and 
acetamide [21j. Figure ^ summarizes the quality of this 
comparison for our final choice of fit parameters, a = 0.6 
and n r = 4.73 x 10" 3 A~ 3 . 



4. RESULTS AND DISCUSSION 
4.1. Pristine surface 

Ref. [T]| reviews in detail the structure of bulk 
Cr203 and the relaxation of its pristine (0001) oxygen- 
terminated surface. Figure |21 shows the relaxed structure 
of our supercell surface slab. The bulk structure consists 
of alternating planar layers of oxygen atoms separated 
by bilayers of chromium. As found in , the primary 
relaxation associated with forming the surface is for the 
oxygen-terminated surface layers to move inward toward 
the bulk crystal with slight in-plane displacements. 

Figure [Ijji.l shows the filled and empty energy lev- 
els from our supercell calculation of the (0001) oxygen- 
terminated surface of Cr2 03 in a vacuum environment. 
Following standard practice, we choose the zero of en- 
ergy to be the Fermi level, the energy below which states 
are fully occupied and above which they are empty. The 
states at the zero line in the figure are thus the high- 
est occupied molecular orbitals (HOMOs) and the first 
states above the line are the lowest unoccupied molecu- 
lar orbitals (LUMOs). The figure shows the levels of the 
pristine surface to be fully filled up to a gap of about 
0.5 eV separating the HOMOs and LUMOs. 




Figure 2. Relaxed structure of pristine surface slab, 
side view ([0001] direction runs vertically up the page): 
oxygen (light grey spheres), chromium (dark grey 
spheres). 
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Figure 3. Near-gap energy levels for (a) pristine, (b) 
hydrogenated and (c) chlorinated system in (1) vacuum 
and (2) dielectric: filled levels (solid circles), empty 
levels (open circles), partially filled levels (crosses). 

To provide a more global view, Figure 21 presents the 
density of states, the number of levels from Figure EJi.l 
per unit energy as a function of energy, computed using 
a Gaussian broadening of width <r=0.41 eV. To under- 
score the distinction between occupied and unoccupied 
states, the figure gives a separate curve for each. Fi- 
nally, as a guide to identification of the bands, the figure 
also contains markers for the LSDA atomic eigenvalues 
of oxygen and chromium, which have been shifted uni- 
formly upward by 4.4 eV to approximately counteract 
the shifting of the Fermi level of the supercell states to 
zero energy. 

Three features in the supercell density of states play 
important roles in the chemistry of this surface. First, the 
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Figure 4. Total density of states of pristine supercell 
slab in vacuum (solid curve) and in solution (virtually 
indistinguishable dashed curve): occupied states (curve 

on left), unoccupied states (curve on right), LSDA 
atomic eigenvalues (vertical arrows) labeled according 
to alignment (|) or anti-alignment (!) with direction of 
net atomic spin. 



highest occupied oxygen orbitals appear as the shoulder 
(from ~-4 eV to ~-2 eV) of the oxygen 2p band, which 
consists of "minority" spin electrons, electrons with spin 
opposite to the net atomic spin. Figure |5ji presents the 
sum of the probabilities associated with states in this 
range and clearly indicates this shoulder to be a surface 
oxygen band. The localization of different spin types to 
either surface corresponds to the fact that the oxygen 
atoms on each surface have opposite net spin direction. 
Next, the highest occupied molecular orbitals (HOMOs) 
overall appear as an sd-hybrid chromium band (from ~- 
2 eV to ~0 eV) consisting of electrons of "majority" spin, 
spin aligned with the net atomic spin. Figure shows 
these states to be bulk chromium states, with atoms al- 
ternating in majority spin direction corresponding to the 
anti-ferromagnetic nature of the bulk material. Finally, 
the LUMOs of the system appear as the low energy ma- 
jority spin shoulder (from ^0 H to ~1 eV) of an empty 
chromium band. Figure shows this shoulder to consist 
of surface states of primarily chromium character pro- 
tected under the outer oxygen layer. This character of the 
HOMOs and LUMOs suggests that any oxygen missing 
from the outer surface would expose a reactive chromium 
layer underneath. 

Upon repeating the pristine surface calculation in the 
presence of a dielectric environment, we find virtually no 
change. There is no change in Figure |21 and, although 
some very small changes are evident in going from Fig- 
ure EJi.l to Figure 01. 2, the global picture of the density 
of states in Figure 0] is visually indistinguishable for vac- 
uum (solid curve) and dielectric (dashed curve). Finally, 
inspection of density maps corresponding to Figures EH3 



t 


99* 


1 








I * 


h t * 








9 i 








* 9 


990 





Figure 5. Contour levels of sum of probability 
densities from the highest energy shoulder of the oxygen 
2p band, side view ([0001] direction runs vertically up 
the page): up spin (black surface), down spin (white 
surface) . 



again shows no noticeable changes. 



4.2. Interaction with hydrogen 

Anticipating bonding with oxygen, we initially placed 
a hydrogen atom in vacuum directly on top of one of 
the surface oxygen atoms, all of which are equivalent by 
symmetry. Figure |BK shows that, upon relaxation, the 
hydrogen atom cants away from the surface perpendicu- 
lar while appearing to form a bond with the underlying 
oxygen atom: the final relaxed O-H distance is 0.95 A, 
quite close to the experimental O-H separation in H 2 
(0.96 A). Figure shows that the canting of the hydro- 
gen atom is in the same direction as one would expect for 
the 2p orbital of the associated oxygen atom given its in- 
plane displacement. Finally, FigurelTUh. which shows the 
total electron density associated with the chemisorbed H, 
confirms the presence of the bond as a small protrusion 
in the density near the hydrogen atom. 

Figure^. 1 shows the filled and empty energy levels of 
the hydrogenated surface in a vacuum environment. As 
with the pristine surface (Figured) > the energy levels are 
fully filled up to a HOMO-LUMO gap, consistent with 
the observed bonding. To better resolve the bond associ- 
ated with the chemisorbed hydrogen, Figure UTI presents 
the local density of states in the vicinity of the hydrogen 
atom, which we compute in the same way as the total 
density of states of Figure 01 but by now weighing each 
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Figure 8. Relaxed structure of surface slab with 
adsorbed hydrogen in (a) vacuum and (b) dielectric: 
oxygen (light grey spheres), chromium (dark grey 
spheres), hydrogen (black spheres). Same view as 
Figure |21 but showing atoms from a single supercell. 



Figure 6. Contour levels of sum of probability 
densities from the chromium sd f band, side view 
([0001] direction runs vertically up the page): up spin 
(black surface), down spin (white surface). 
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Figure 7. Contour level of sum of probability 
densities from the lowest energy shoulder of the 
chromium sd j spin band, side view ([0001] direction 
runs vertically up the page): up spin (black surface), 
down spin (white surface). 



Figure 9. Relaxed structure of single supercell with 

adsorbed hydrogen in vacuum, top view ([0001] 
direction normal to page): oxygen (light grey spheres), 
chromium (dark grey spheres), hydrogen (black sphere). 



level with the probability of an electron in the level being 
within 0.69 A of the proton. The local density of states 
shows that the hydrogen atom interacts mostly with the 
surface oxygen 2p band. A plot of the total density asso- 
ciated with this surface band, Figure fT5k. confirms that 
it contains most of the density protrusion associated with 
the O-H bond. Finally, FigureEt shows that the HOMO 
of the hydrogenated surface, while maintaining signifi- 
cant bulk chromium character, indeed localizes near the 
hydrogen atom. 

To determine the final atomic configuration in the pres- 
ence of a solvent, we begin with the positions from Fig- 
ure |SJi and relax the atomic coordinates within the ap- 
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(a) 



Figure 10. Contour level of 0.68 e~/A 3 total electron 

density for hydrogen atom adsorbed on 
oxygen-terminated (0001) surface in (a) vacuum and (b) 
dielectric, side view ([0001] direction runs vertically up 
the page). 
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Figure 12. Contour level of sum of probability 
densities associated with oxygen 2p surface band in (a) 
vacuum and (b) dielectric: up spin (black surface), 
down spin (white surface). 
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Figure 11. Local density of states within 0.69 A of 
the proton for hydrogen interacting with 
oxygen-terminated (0001) surface in vacuum (solid 
curve) and solution (dashed curve). Same conventions 
as Figure 0] 



proximate joint density functional (l.'il II) until the maxi- 
mum force on any atom is less than 0.3 eV/A. Figure^, 
which displays the resulting configuration, shows that the 
presence of the dielectric has a dramatic effect on the 
hydrogen. The nearest oxygen-hydrogen distance has in- 
creased to 2.3 A, clearly breaking the O-H bond and re- 
turning the hydrogen atom to the solution. Consistent 
with this picture, the largest residual force remains on 
the hydrogen atom in the direction tending to push it 
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Figure 13. Contour level of HOMO of hydrogen 
interacting with the oxygen-terminated (0001) surface 

in (a) vacuum and (b) dielectric, side view ([0001] 
direction runs vertically up the page): up spin (black 
surface), down spin (white surface). 



yet further from the surface. Lack of any indication of 
the presence of the hydrogen atom in the resulting to- 
tal charge density, Figure 110b . indicates that the atom 
enters the solution as an ion. 

Figure |2t>. 2 shows that upon the removal of the pro- 
ton from the surface, the excess electron from the O-H 
bond appears simply as a donated conduction electron 




Figure 14. Relaxed structure of surface interacting 
with chlorine in vacuum: oxygen (light grey spheres), 
chromium (dark grey spheres), chlorine (black sphere): 
top view ([0001] direction perpendicular to page). 

just above the energy gap. Consistent with this dona- 
tion, the dashed curve in Figure 1111 shows a dramatic 
reduction in the local density of states near the proton, 
and the HOMO in Figure 113b now shows more of the 
surface chromium character of the original LUMO band 
from Figured 

4.3. Interaction with chlorine 

Anticipating the possibilities of ionic bonding for chlo- 
rine, we initially placed (in the lxl surface supercell) a 
chlorine atom in vacuum directly above each of the two 
distinct types of Cr 3+ site from the outermost chromium 
bilayer and found the site above the innermost of the 
two component layers to be favored by 0.4 eV. Figure PHI 
shows a top ([0001]) view of the relaxed configuration 
for this site when computed within the 4x4 supercell. 
We find relatively little relaxation from the clean sur- 
face structure (movement of all surface atoms is less than 
0.01 A) with the chlorine ion settling upon relaxation to a 
position with a chlorine-oxygen separation of 2.6 A, only 
20% smaller than the sum of the nominal ionic radii, 
3.1 A. Fi gure 115b shows a side view of the total elec- 
tron density for a surface with adsorbed CI in vacuum, 
illustrating the physisorbed nature of the interaction. 

Turning to the energy levels, Figure |2t-l shows the 
levels near the gap. In this case, three two-thirds filled 
states (degenerate to within 19 meV~ 220 K) appear just 
below the top of the gap, indicating the presence of a hole 
in the Cr band, which we interpret as arising from the 
chlorine atom absorbing an electron from the chromium 
oxide to become Cl~. To explore local effects from the 
adsorbed chlorine, Figure ED presents the local density 
of states, weighing each state with the probability of an 
electron being within 1.6 A of the chlorine nucleus. In 
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Figure 15. Contour level of 0.68 e~/A 3 total electron 

density for chlorine atom adsorbed on 
oxygen-terminated (0001) surface in (a) vacuum and (b) 
dielectric, side view ([0001] direction runs vertically up 
the page). 

contrast to the local density of states of the hydrogen 
calculation, the appearance of the oxygen band is signif- 
icantly reduced and there is much stronger mixing with 
the bulk chromium states. This mixing corresponds to 
alignment of the barely bound Cl~ states with the bulk 
HOMO chromium band as the chlorine ion draws elec- 
trons from the bulk chromium band, which is serving as a 
reservoir of electrons. Finally, Figure ITVk shows the sum 
of the electron probabilities in the partially filled states 
at the Fermi level, which are thus both the HOMOs and 
the LUMOs. Interpreted as the sum states lacking ex- 
actly one electron from full occupancy, the figure shows 
the spatial distribution of the hole which the formation 
of the Cl~ injects into the chromium-oxide slab. As one 
would expect, this (positive) hole tends to localize to the 
vicinity of the Cl _ ions. 

Upon relaxation of the physisorbed chlorine surface in 
the presence of the solvent, we find there to be very lit- 
tle relaxation (no more than 0.01 A for any atom), little 
difference in the total charge density fFigure ITS'b'). no 
change in the presence of partially filled states at the 
Fermi level (Figure ISfc. 2), and very little difference in the 
local density of states (Figure ^jj) or the spatial distribu- 
tion of the hole injected into surface, Figure IT7b. 

4.4. Conclusions 

Above, we introduce the novel approach of using a joint 
density-functional theory to treat an ab initio electronic 
structure calculation in the presence of a liquid solvent 
such as water. The resulting approach is the first practi- 
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Figure 16. Local density of states within 1.6 A of the 
chlorine nucleus for chlorine interacting with 
oxygen-terminated (0001) surface in vacuum (solid 
curve) and solution (dashed curve). Same conventions 
as Figure 0] 
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Figure 17. Contour level of HOMOs and LUMOs 
(equivalent in this case) of chlorine interaction with 
oxygen-terminated (0001) surface in (a) vacuum and (b) 
dielectric, side view ([0001] directions runs vertically up 
the page). 



cal ab initio approach for the treatment of the interaction 
of complex surfaces with chemical species in the presence 
of a dielectric environment. 

Through this approach, we find the mode of interac- 
tion of the oxygen-terminated CT2O3 (0001) surface with 
hydrogen to be covalent bonding while that with chlorine 
to be ionic bonding. The presence of a dielectric solvent 
has very little effect on the pristine surface or on its inter- 



action with chlorine, while it has a dramatic affect on the 
interaction with hydrogen. In vacuum, hydrogen readily 
forms an O-H bond with the outermost layer of atoms of 
the surface. In the presence of water, the strong screen- 
ing associated with dielectric effects in the vicinity of the 
proton (ultimately via hydrogen bonding interactions) so 
weakens the attractive potential of the proton that the 
covalent bond is broken, the electron is released into the 
surface and the proton solvates. 

In contrast, the interaction with chlorine in vacuum is 
already ionic, with a neutral chlorine atom having suffi- 
cient electronegativity to draw an electron from the bulk 
of the crystal and thus ionize to CI - while injecting a 
hole into the bulk. The presence of a dielectric solvent 
tends to screen the excess charge on the chlorine ion, 
thereby only further stabilizing this form of interaction 
with the surface so that there is little change in this case 
when going from vacuum to a dielectric environment. It 
thus appears that the primary reason why the solvent has 
a much greater effect on the interaction with hydrogen 
than with chlorine is that a dielectric environment gener- 
ally favors the formation of ions and the surface interac- 
tion with chlorine is already ionic whereas the interaction 
with hydrogen in vacuum is covalent. 

We believe that there would be little effect on these 
conclusions were the calculations to be performed with 
ions rather than atoms. Doing so would involve prin- 
cipally removing a single electron from the calculation 
for each hydrogen atom or adding an electron for each 
chlorine atom. For the chlorine cases, this would sim- 
ply remove the relatively delocalized hole from the bulk 
chromium band and thus likely have little effect on the 
final results. For the interaction with hydrogen, the re- 
moval of an electron would, in the dielectric case, simply 
remove the relatively delocalized donated electron or, in 
the vacuum case, likely simply introduce a relatively de- 
localized hole into the chromium band. In either of these 
cases for hydrogen, we again would expect little disrup- 
tion of the chemical integrity of the surface. 

Overall, a novel picture emerges to explain how the 
oxygen-terminated surface of Cr203 is particularly pro- 
tective in a hydrochloric acid solution. The outer oxy- 
gen layer provides a natural barrier to interaction with 
chlorine but might be expected to interact strongly with 
protons. However, through dielectric screening effects, 
it is the aqueous environment itself which eliminates the 
outer oxygen layer's natural tendency to interact with 
hydrogen. 

These first calculations of surface chemistry in the 
presence of a solvent make clear the need for additional 
work to complete the picture of the passivating effects 
of chromium oxide. We would next like to study the 
interaction of chlorine and hydrogen with a chromium- 
terminated surface, whose HOMOs would then be ex- 
posed on the surface rather than protected under the 
outer oxygen layer. We also would like to explore pit 
corrosion by calculating the interaction of the aforemen- 
tioned species with a step in the passivating oxygen layer. 
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Finally, we would like to contrast the interactions of chlo- 
rine in such systems to the interactions of fluorine and 
bromine in order to better understand the corrosive suc- 
cess of chlorine relative to these other species. 
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